Extended Poisson— Boltzmann descriptions of the electrostatic double layer: 
implications for charged particles at interfaces 
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I. INTRODUCTION 

Whenever a charged surface is immersed in a thermal 
solvent with mobile co- and counterions (usually an aque- 
ous electrolyte), a screening layer of counterions will form 
in the vicinity of the charged surface, giving rise to the 
so called "double layer" . Systematic investigations of the 
double layer physics and chemistry have been conducted 
for more than one hundred years and have been moti- 
vated by various practical questions such as the behavior 
of battery electrodes or the stability of charged colloid so- 
lutions. Very often surfaces can be highly charged (e.g. 
colloid surfaces attain charge densities of about one el- 
ementary charge e per nm 2 ) such that the electrostatic 
potential ip near the surface exceeds the equivalent of the 
thermal energy ksT = 1/(3 by far (for room temperature 
k^T/e corresponds to a potential of 25 mV). Thus the 
arrangements of water molecules and counterions near 
and at the surface should display fairly strong correla- 
tions, resulting in e.g. nontrivial potential and dielectric 
profiles. 

Such correlations are neglected in the "work horse" 
model of charges in solutions, the Poisson-Boltzmann 
theory. Here, Poisson's equation V(eVV'(r)) = 
— e(c+(r) — c_(r)) is coupled with a strongly simplify- 
ing assumption on the density distribution of positive 
and negative ions, c± — c s exp(=pe/?V'(r)), i.c a Boltz- 
mann distribution of charged point particles (having no 
correlations) subject to the electrostatic potential (the 
ions are monovalent for simplicity and e is the dielectric 
constant of the solvent, c s is the bulk number density 
of positive/negative ions in the solvent). The resulting 
Poisson-Boltzmann equation for the electrostatic poten- 
tial 



V • (eVV>) = 2c s esinh(/3eV>) 



(1) 



still poses a formidable nonlinear problem, analytically 
solvable only in special circumstances such as the ge- 
ometry of a charged wall. For e = const., it is conve- 
nient to introduce the Debye-Hiickel screening length via 
= (e/(2c s /3e 2 )) 1//2 and the dimensionless potential 
<t> = fieij). Using these quantities, the Poisson-Boltzmann 
equation becomes V • V</> = « 2 sinh0, and the signifi- 
cance of the screening length becomes immediately ob- 



entails an exponentially decaying potential with charac- 
teristic length ft -1 around a localized charge distribution. 

Very often charged particles in electrolytes are char- 
acterized using properties of the linearized Poisson- 
Boltzmann equation. The effective or renormalized 
charge or charge density is the one appropriate for 
the charged particle immersed in a hypothetical "linear 
Poisson-Boltzmann" medium which exhibits the same 
asypmptotic potential decay as the "real" charged par- 
ticle in the "real" (and complicated) electrolyte. If two 
particles immersed in the electrolyte are sufficiently far 
apart, then their interaction energy is proportional to 
the product of their effective charges - this illustrates 
that the charge renormalization concept is very useful in 
characterizing charged colloid suspensions. However, the 
effective charge density of particles in bulk electrolytes 
contains very little information about the precise struc- 
ture and potential profile in the double layer. This is 
best illustrated by the well-known solution for a charged 
wall immersed in an electrolyte in Poisson-Boltzmann 
theory. Using dimensionless charge densities given by 
a = (3e/(Ke) a the effective charge density a of the 
wall is given in terms of the bare charge density a c by 



a cff = 4 



1 + ^ 



(2) 



vious in its linearized form, (V ■ V 



0, which 



For colloidal matter, <r c < 1 e/nm 2 and except for situ- 
ations where k _1 < 0.3 nm (or c s > 1 M) we are in the 
saturated regime where a cS = 4. We can draw an im- 
portant conclusion from this result: Although we might 
expect the Poisson-Boltzmann description to fail near 
the charged particles for typical colloid charge densities, 
we must expect that the dimensionless effective charge 
density is roughly constant for modest variations of a c 
and k _1 and its determination will not allow to discrim- 
inate between different models of the double layer. Thus 
by measuring pair interactions between charged particles 
we mostly learn which area A of the particle's surface is 
covered with charge (since q cS = Aa cS ) and we do not 
have access to, say, the potential at the particle surface. 
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II. CHARGED PARTICLES AT ELECTROLYTE 
INTERFACES 




FIG. 1. Sketch of electric field lines around a charged col- 
loid located at the interface between water and a nonpolar 
medium. The asymptotics of the electrostatic potential and 
the associated electric field at the interface is mainly deter- 
mined by field lines originating from the colloid charges and 
passing through the colloid and the nonpolar medium. Field 
lines through the water (containing ions) are screened (indi- 
cated by the dashing of the lines). 

The situation is different for charged particles at the 
interface between an electrolyte and a nonpolar medium 
(air, oil), see Fig. [TJ If the dielectric constant of the 
nonpolar medium is low (e/eo < 2, e is the dielectric 
constant of vacuum) then surface charges can only reside 
on the colloid surface exposed to the electrolyte. The 
double layer screens the colloid charge only to the ex- 
tent that no monopole field is present along the inter- 
face. However, higher multipole fields are present, and 
the leading multipole is a dipole with moment p z oriented 
perpendicular to the interface (regardless of the charge 
distribution on the colloid [H). Consequently, the elec- 
trostatic interaction energy between two charged colloids 
at the interface with mutual distance d is long-ranged 
and ocpl/d 3 . If the local curvature at the colloid surface 
is smaller than k then the double layer can be approxi- 
mated by the double layer of a charged wall in the vicin- 
ity of each point at the colloid-electrolyte interface. The 
dipole moment p z = A cS p' wall can be calculated from the 
aerial dipole moment density p' wal i and an effective area 
A eS which depends on the precise colloid geometry. The 
crucial observation is that p' u = eoVwu is proportional 
to the contact potential at the wall.[? ] Thus the inter- 
action energy at large distances between two colloids is 
proportional to the square of the wall contact potential, 
and such a dependence is very much different from the 
dependence on the effective (saturated) charges for the 
interaction between colloids in bulk electrolyte @. The 
Poisson-Boltzmann equation gives j3ep' waa « 2eo In a c for 
a c > 1, so already here is seen that the dipole moment 
will not saturate with increasing charge density on the 
colloid and depends logarithmically on both screening 
length and charge density, p z ex m(K _1 cr c ) + const, (see 
the definition of a). 

Experiments. - The dipole moment can be calculated 
from the interaction potential between charged colloids 



and this can be inferred through laser tweezers on isolated 
pairs of colloids, from inversion of measured pair corre- 
lations or from elasticity measurements on 2d crystals at 
the interface. Whereas the logarithmic dependence on 
k _1 has received much support from laser tweezer data 
pl| . the absolute magnitude of p z from the Poisson- 
Boltzmann equation is too small by about a factor of 
5 compared to results of careful measurements using all 
three routes [llj . This points to a severe underestimation 
of the contact potential at a charged wall in the Poisson- 
Boltzmann equation. 

III. MODIFICATION STRATEGIES FOR 
IMPROVING THE POISSON-BOLTZMANN 
DESCRIPTION 

The deficiencies of and possible amendments to the 
Poisson-Boltzmann equation ([T]) have been discussed for 
a long time, starting with the work of Bikerman 0). 
Recently, starting with Ref. Q, the modified Poisson- 
Boltzmann approach received renewed interest. We are 
interested in modifications which preserve the type of 
equation as a differential equation for the potential if. 
The electrostatic interactions thus remain on the mean- 
field level, and the modifications introduce contributions 
that do not stem from correlations. This approach is 
suitable to investigate certain aspects of "ion specificity" 
in the weak-coupling limit. 

1. Finite size effects. - A finite size of the co- and 
counterions (e.g. through a hard-sphere volume 
v) will inevitably lead to layering effects in the 
density profile of the ions near the charged sur- 
face. In aqueous solution, the effective ion size 
is increased as a result of hydration. Thus, in 
a strict sense, the finite ion size contributions re- 
quire nonlocal treatment, however, one can cap- 
ture these effects using local corrections, for exam- 
ple, by taking into account a local change in the 
ion chemical potential [i± (or ion osmotic pressure 
p) through local density changes. In the Poisson- 
Boltzmann equation ions are treated as ideal, thus 
/3fj,± = lnvc± (or /3p = p). A popular modification 
(MPB, modified Poisson-Boltzmann equation) Q 
has been derived by lattice gas considerations, lead- 
ing to /3/x± = \avc± — ln(l — vc+ — vc-) (or /3p = 
— v^ 1 ln(l — vc+ — vc-)), where the excess chemical 
potential, v -1 ln(l — vc + — wc_), accounts for the re- 
pulsive (excluded volume) ionic interactions. Such 
a local correction places an upper bound on the lo- 
cal density, leading to density "saturation" in the 
region of high ionic concentration. Far from the 
concentrated region, the density is merely shifted 
by some length d, c MPB (x) = c PB (x — d), in re- 
lation to the Poisson-Boltzmann result. Thus, the 
screening behavior far from the concentrated region 
is not exactly altered, but rather the effective posi- 
tion of a charged surface is shifted into the solution. 
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Roughly, the changes associated with the finite size 
effects can be quantified with an additional length 
scale vcr c , which is the size of a layer of counterions 
that would form if all counterions were forced into 
close-packed configuration. The extent of the finite 
size effects is estimated as the ratio of this length 
with other relevant lengths, such as the screening 
or Gouy-Chapman length. Although the lattice 
gas modification corresponds to a not very precise 
equation of state for hard spheres, application of 
MPB to inhomogeneous ionic solution yields pre- 
dictions that compare rather well with exact re- 
sults. This points to some internal cancellation of 
errors. 

2. Dielectric inhomogeneity of the solvent. - The sol- 
vent medium, made up of water molecules with 
permanent dipoles, acts as a dielectric which is 
polarized in an external electrostatic field and 
thus screens the electrostatic potential within the 
medium. Thermal fluctuations favor isotropic ori- 
entation of dipoles, and the electrostatic interac- 
tions favor alignment with the field lines. As long as 
the thermal energy exceeds the electrostatic coun- 
terpart, the polarization is proportional to an exter- 
nal field and the dielectric constant is field indepen- 
dent. But if the two energy scales become compara- 
ble, polarization saturation takes place: the dipoles 
reach perfect alignment with the field lines and no 
longer dissipate electrostatic energy above a cer- 
tain threshold. This behavior is described through 
the Langcvin function C(x) = cotha; — 1/x, and 
the polarization with nonlinear dependence on the 
field strength is P = p Cd£{(3po\ V-01), where p is 
the dipole moment of a solvent molecule and Cd is 
the dipole density. This amounts to a change in 
dielectric constant according to 

£M) 
e e ff -> eo + PoCd r=-rj ■ (3) 

Here, it is assumed that the solvent is incompress- 
ible and the modification leads to the Langevin 
Poisson-Boltzmann equation (LPB) Q- For small 
electrostatic field, e ~ e + /3pg c d/3 and in the limit 
of large field e — > eo . Uniform external fields affect 
only the dipole orientation. Nonuniform fields ex- 
ert, in addition, a translational (dielectrophorctic) 
force, so that a dipole migrates toward regions with 
stronger fields. If one treats the solvent as consist- 
ing of point dipoles, as in Ref. [J (DPB, dipolar 
Poisson-Boltzmann equation), this force leads to 
the strong dependence of the solvent density on 
the field, Cd{fipo\ V"0|), with the functional form 
Cd{x) = smhx/x. However, due to hydrogen bonds 
which lead to a high density of water at ambient 
conditions, the assumption of incomprcssibility is 
more realistic and the resulting e qualitatively re- 



produces the drop with increasing field strength as 
seen in the simulations [l2j ■ 

3. Ionic polarizabilities. - Ions can further be charac- 
terized by including their polarizability a, so that 
they acquire an induced dipole moment in an ex- 
ternal field. Part of the electrostatic energy is dis- 
sipated when it polarizes ions. This leads to a di- 
electric constant with a term linear in the ion con- 
centration, 

e cff e + a(c + +c_). (4) 

The presence of induced dipoles, in turn, con- 
tributes to a dielectrophorctic interaction of an ion 
with a field. Accordingly, the density distribution 
of ions becomes [j| , 

c± =c s cxp(Te/# + /fo|V^| 2 /2). (5) 

The above theoretical framework can be used to de- 
scribe a somewhat different but related situation. 
A dissolved charged particle (for simplicity taken to 
be a point) is tightly bound to solvent dipoles form- 
ing its hydration shell so that the dipoles become 
unresponsive to presumably a weaker external field. 
This produces a dielectric cavity which reduces the 
dielectric constant of a solvent, on the one hand, 
and leads to a migration of a hydrated ion away 
from a strong field, on the other hand. These trends 
are captured with polarizable model where the "ef- 
fective" polarizability is a negative valued quantity 
0). Interestingly enough, both negative and posi- 
tive a lead to a spreading out of a diffuse double 
layer. 

4. Stern layer. - A layer of solvent in direct contact 
with an interface has different properties from the 
bulk solvent. First of all, ions cannot penetrate 
into it due to their finite size. Second of all, the 
solvent within this layer tends to be more struc- 
tured so that its dielectric response is expected to 
be suppressed in relation to that in the bulk. A fre- 
quent solution to incorporating this contribution is 
to introduce the Stern layer adjacent to a sharp in- 
terface with the width corresponding to a radius of 
a bare or hydrated ion and the dielectric constant 
the same as or lower than the bulk value. The in- 
troduction of this layer does not change anything 
far from an interface, but it increases the contact 
potential, ip c — > ip c + a c R/ e s tcrn, where R is the 
radius of an ion and e s tem is the dielectric constant 
with the Stern layer Q . Incorporation of the Stern 
layer at an air-water interface was found to give an 
accurate predictions for the excess surface tension 
of strong electrolytes (l3j . 

As an example, in Fig. [5] we present a comparison be- 
reen PB and MPB, LPB (routes (1) and (2) from above) 
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for the contact potential dependence at a charged wall on 
the charge density. Finite ion size and dielectric satura- 
tion lead to a noticeable increase in the contact potential 
which is in line with the experimental results (discussed 
above) for the dipole moment of charged colloids at in- 
terfaces. 

Fig. [3] compares countcrion density profiles for various 
models. The "saturation" effect is clearly seen for the 
MPB model. The mechanism behind the spreading-out 
of the diffuse layer for the PPB model is connected with 
the dielectrophoretic repulsive force of ions with negative 
polarizability. Only the LPB model causes a stronger at- 
traction toward the wall. This is caused by the lowering 
of the dielectric constant close to a wall. The combina- 
tion of effects would require an appropriately combined 
model, an example is the combination of LPB and MPB 
as seen in the thick green dot-dashed curve in Fig. [5] 
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FIG. 2. Contact potential (proportional to a colloid dipole 
moment at an interface) for PB and LPB, with and without 
steric corrections due to finite ion size (v — a 3 ). The symbols 
correspond to numerical data points and dashed lines corre- 
spond to analytical expressions derived in Ref. Q. Clearly, 
the contact potential in PB is noticeably smaller than in the 
modified theories for higher charge densities, thus dipole mo- 
ments will be underestimated by PB alone. 



IV. OUTLOOK TO NONLOCAL 
DESCRIPTIONS 

In the models outlined above, the electrostatic inter- 
actions arc described on the mean-field level, and the 
repulsive-core interactions among hydrated ions and sol- 
vent molecules arc incorporated as local contributions. 
These models, even if quantitatively not precise, are able 
to capture correct trends. The drawback of the local ap- 
proximation for the excluded volume interactions is that 
the "saturated" density profile is an unphysical artifact 
of the method itself. This leads to a violation of the 
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FIG. 3. Counterion density near a charged wall for the mod- 
els discussed in Sec. IIIII The common parameters are: a c = 
0.4 C/m 2 , c s —0.1 M, and As = 0.7 nm. Specific parameters 
for each model are: for MPB the ion radius is 0.35nm, for 
LPB (point ions) the dipole concentration is d = 55M (as 
in water) and the permanent dipole moment is po = 4.85 D, 
for PPB the polarizability is a/(47re ) = -8 M -1 . Only LPB 
leads to an increased concentration near a wall, while the 
other contributions tend to spread out a diffuse layer. 



contact value theorem (the relation between the density 
at the wall, solvent pressure and wall potential). Despite 
this, the electrostatic properties at a charged wall and far 
from the wall are captured with reasonable accuracy, de- 
spite the unphysical details in-between. A possibly more 
accurate treatment of the repulsive-core interactions can 
be achieved via density functional theory. In the example 
of Ref. , ions are treated as charged hard spheres using 
the accurate fundamental measure theory, but for sim- 
plicity in the background of a structureless solvent. This 
corresponds to a nonlocal extension of the MPB modifica- 
tion discussed above. Finally, the mean-field treatment 
of electrostatics neglects correlated fluctuations, which 
not only correct, but, in the so-called strong-coupling 
limit, give rise to altogether new phenomena, such as 
the attraction between the same-charged macroparticles 
[Til . [l5| . While the applicability of the mean- field treat- 
ment breaks down in the strong-coupling limit, the cor- 
rect theoretical treatment of correlations continues to be 
a challenge. This challenge is increased if both the corre- 
lations due to the short-range solvent and the long-range 
electrostatic interactions are important. Currently, the 
solutions are looked for using either the density func- 
tional theory approach or the field-theoretical methods. 
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